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I review some recent progress, leading to a substantial reduction in the number of non- 
perturbative parameters, in the calculation of inclusive quarkonium decay widths in the frame- 
work of non-relativistic effective field theories. 

<n : 

£>V 1 Introduction 

& : 

^ | Considered as bound states, heavy quarkonia (e.g. T, Xb, Xc •••) have two main character- 

^r) ' istics: 1) a mass scale, which is in the perturbative regime: mj ~ 5 GeV, m c ~ 1.5 GeV; 2) 

level splittings of a relative size that is typical of non-relativistic systems. Being non-relativistic 
bound states, heavy quarkonia are characterized by at least three hierarchically ordered energy 
J> ' scales: the hard scale m, the momentum scale mv and the bound-state energy scale mv , where 

v <C 1 is the heavy-quark velocity inside the bound state. H 

The different scales entering the quarkonium dynamics may be systematically integrated 
\ out, leading from QCD to simpler but equivalent effective field theories (EFTs). NRQCD is the 

EFT obtained by integrating out the hard scale m. H'i This m being larger than the scale of non- 
perturbative physics, Aqcd, the matching to NRQCD can be done order by order in a s . Hence, 

^^^^ ' " (dri) 

the NRQCD Lagrangian can be written as a sum of terms like f n O n n /m d ~ 4 , ordered in powers 
Q_i| of a s and v. More specifically, the Wilson coefficients f n are series in a s (m) and encode the 

Q_i| ultraviolet physics that has been integrated out from QCD. The operators O n of dimension d n 

describe the low-energy dynamics and are counted in powers of v. Heavy quarkonium decays are 
processes that take place at an energy-transfer scale of the order of the heavy-quark mass. Their 
. ^ \ signature is encoded in the imaginary part of the Wilson coefficients of the 4-fermion operators 

(O n dn) = ^K n xx ] K' n ^) in the NRQCD Lagrangian. The NRQCD factorization formula for 
■ quarkonium inclusive decay widths into light hadrons (LH) readsB 

2Tm f 

T(H - LH) = Y, (H|^„XX^|H). (1) 

The 4-fermion operators are classified with respect to their rotational and spin symmetry (e.g. 
0( +1 Sj), 0( 2S+1 Pj), ...) and of their colour content (octet, Og, and singlet, Oi, operators). 
Singlet operator expectation values may be easily related to the square of the quarkonium wave 
functions (or derivatives of it) at the origin. These are unknown non-perturbative parameters 
and may be fitted from experimental data, calculated on the lattice or taken from phenomeno- 
logical potential models. Octet operator expectation values are also unknown non-perturbative 
parameters to be fitted from the experimental data, or to be calculated on the lattice. 

NRQCD still contains dynamical degrees of freedom associated with energy scales larger 
than the ultrasoft scale mv 2 . @ j] Hence, pushing further the EFT programme for non-relativistic 
bound states, further simplifications occur if we integrate out those degrees of freedom. We call 
pNRQCD the resulting EFT. We will consider pNRQCD under the condition Aqcd S> mv 2 . 
Then, two situations are possible. First, the situation when mv ^> Aqcd ^ mv 2 . In this case 

a As a consequence the power counting of the NRQCD operators is not unique, since they depend, in general, 
on the scales mv, mv 2 and Aqcd- 



the soft scale mv can be integrated out perturbatively. This leads to an intermediate EFT 
that contains singlet and octet quarkonium fields and ultrasoft gluons as dynamical degrees of 
freedom. The octet quarkonium field and the ultrasoft gluons are eventually integrated out by 
the (non-perturbative) matching to pNRQCD.i Second, the situation when Aqcd ~ mv ti In this 
case the (non-perturbative) matching to pNRQCD has to be done in one single step. □ Under 
the circumstances that other degrees of freedom develop a mass gap of order Aqcd ° r that they 
play a negligible role, the quarkonium singlet field S remains as the only dynamical degree of 
freedom in the pNRQCD Lagrangian, which readsi'ii £ p nrqcd = Tr (ido — 7i)S\, TC being 
the pNRQCD Hamiltonian, to be determined by matching pNRQCD to NRQCD. The inclusive 
quarkonium decay width into light hadrons is given by 

r(H — > LH) = —21m (n, L, S, J\H\n, L, S, J), (2) 

where |n, L, S, J) is an eigenstate of Tt with the quantum numbers of the quarkonium state H. 
From the matching we obtain schematically: 

= E + {^(r), A} £ ^B n + V^(r)V £ + . . . , (3) 

n n n 

where f n are the matching coefficients inherited from NRQCD, and A n , B n , ... are non- 
perturbative operators, which are universal in the sense that they do not depend either on the 
heavy-quark flavour or on the specific quantum numbers of the considered heavy-quarkonium 
state. Inserting Eq. (|3|) into (||) and comparing with Eq. (|l|), we see that all NRQCD matrix 
elements, including the octet ones, can be expressed through pNRQCD as products of the uni- 
versal non-perturbative factors by the squares of the quarkonium wave functions (or derivatives 
of it) at the origin. This drastically reduces the number of unknown non-perturbative factors 
to be introduced in order to describe the whole set of charmonium and bottomonium inclusive 
decays and makes new theoretical predictions possible. As an example, in the following we will 
discuss the case of P-wave inclusive quarkonium decays into light hadrons. Q 



2 P-wave decays in NRQCD and pNRQCD 



In NRQCD the P-wave inclusive decay width for the I 
states is given at leading (non-vanishing) order in v by: 



(h) and S = 1 (x) quarkonium 
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(4) 



1 Rp H n \h \0 8 { Sq) \h), tor J = 0,1, 2 (5) 

where R' P is the derivative of the P-wave function at the origin. We stress that, according to 
the power counting of NRQCD, the octet contribution (/i|Og( 1 S'o)|/i) is as relevant as the singlet 
contribution; the above formula may by no means be used without considering it (for instance 
in extracting a s (m), as can still be found in some recent literature). We also recall that the 
octet contribution reabsorbs the dependence on the infrared cut-off \i of the Wilson coefficients 
Im/i(P). Finally, we see that in NRQCD the 8 P-wave bottomonium states (IP, 2P) and 
the 4 P-wave charmonium states (IP), which lie under threshold, depend at leading order in 
the velocity expansion on 6 non-perturbative parameters (3 wave functions + 3 octet matrix 
elements). 

In pNRQCD the P-wave inclusive decay widths are given at leading order in v by:! 
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Figure 1: Plot of the 1-loop RG-improved expression of £ vs. fj,: £(fj,) 



£(m) + 7p In s 7 ) . £(m) has been 



extracted from the charmonium P-wave data. The error band accounts only for the uncertainties inherited from 

the charmonium data. 
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for J = 0,1, 2 



(7) 



where £ = -j dt t 3 {gE a (t, 0)$ a6 (i, 0; 0)fifE b (0, 0)^ is the universal non-perturbative param- 
eter that describes P-wave quarkonium decays in pNRQCD. A plot of £ as a function of the 
factorization scale \i is shown in Fig. |l[ In order to better understand the meaning of £, we can 



express it in terms of the eigenstates \n) and the energies E n of the static QCD Hamiltonian: 

(<%E|n) • (n| 5 E|0) 



£ = 18 
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bmce £/q ~ 

excitations we have assumed a mass gap of order Aqcd), we see from Eq. (H) that the correlator 
£ resums all the corrections (mv 2 / A-qcd) u , which one would expect to show up in pNRQCD 
after having integrated out the scale Aqcd- 

By comparing Eqs. (Q) and (|B|) with Eqs. @ and (0) we get at leading order in v the 



(it is the quarkonium static energy) and E^1 Q ~ Aqcd (for higher gluonic 
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relation between the octet matrix element of NRQCD and £: (/i[Og( 1 S'o)|/i) 
The quarkonium-state dependence factorizes in the pNRQCD formulas. This allows some new 
predictions with respect to NRQCD, which are synthetized by the formula (valid at leading 
order in v): 

r(H( 2S+1 nPj) -> LH) 811m/ 1 ( 25+1 Pj) +Im/ 8 ( 2S+1 S 



\£ 



r(H( 25 ' +1 nPj/) — > LH) 8Umf 1 ( 2S '+ 1 Pj,)+Imf 8 { 2S, + 1 Ss>)£ : 



(9) 



where the left-hand side is a ratio between inclusive decay widths of P-wave quarkonia with the 
same principal quantum number n and the right-hand side no longer depends on n and has the 
whole flavour dependence encoded in the Wilson coefficients, which are known quantities. 

In practice, the 12 P-wave quarkonium states, which lie under threshold, depend only, 
in pNRQCD at leading (non-vanishing) order in the velocity expansion, on 4 non-perturbative 
parameters (3 wave functions + 1 chromoelectric correlator £). The reduction by 2 in the number 
of unknown non-perturbative parameters with respect to NRQCD, allows us to formulate two 
new statements. Using the values of £ from Fig. [l] we get (at NLO): 

r( X 6o(lP) -> LH) _ r( X60 (2P) - LH) 



r(x W (lP) - LH) r(xti(2P) - LH) 



.0± 1.3, 



(10) 
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Figure 2: The left-hand side of Eqs. (JToJ) and ([H]) plotted vs. fi. We have used Eq. (^) and the values of £ from 
Fig. |lj The LO and NLO bands refer respectively to the Wilson coefficients at leading and next-to-leading order. 



or alternatively 

r( X6 i(lP) - LH) r( X6 i(2P) - LH) 50 +o.06 ( , , , 



r(x6 2 (iP) -> LH) r( X62 (2P) -> LH) 
The errors refer only to the uncertainties displayed in Fig. [I]. In Fig. |2| we plot the above ratios 
of decay widths as functions of the factorization scale fi. The figures show a stable result. 



3 Conclusions 

By integrating out degrees of freedom associated with energy scales larger than mv 2 , octet matrix 
elements in NRQCD can be written as products of wave functions at the origin by universal 
non-perturbative factors. This reduces the number of unknown matrix elements of NRQCD and 
enables definite predictions for inclusive decay ratios of states with different flavour or principal 
quantum number. As an explicit example we have shown the case of P-wave quarkonium 
inclusive decay into light hadrons. The same program may be carried out for 5-wave decays, i 
At 0(mv 5 ) the hadronic inclusive decay widths of the 10 S-wave quarkonium states, which lie 
under threshold, are described in NRQCD by about 30 non-perturbative matrix elements: 10 
wave functions and 20 further matrix elements. These 20 matrix elements reduce dramatically 
in pNRQCD to few universal correlators, which may be fixed phenomenologically or by lattice 
calculations. Noting that in the hadronic/electromagnetic ratios the wave- function dependence 
drops out, we expect that in the near future pNRQCD will enable QCD to access a large number 
of new results concerning heavy quarkonium inclusive and electromagnetic decays. 
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